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Abstract 

This article introduces the subject of quasi-local horizons at a level suitable for physics 
graduate students who have taken a first course on general relativity. It reviews properties of 
trapped surfaces and trapped regions in some simple examples, general properties of trapped 
surfaces including their stability properties, the definitions and some applications of dynamical- 
, trapping-, and isolated-horizons. 

1 Introduction 

The first conception of a black hole was due to Michell and Laplace in the 18**" century. They 
viewed it as a star whose gravitational field is so strong that the Newtonian escape velocity 
y/2GMjR (with M and R being the mass and radius of the star respectively) is larger than 
the speed of light. The condition on the escape velocity leads to the inequality R < 2GM/c^ 
which, remarkably, also holds in general relativity. While such a star would have the property 
that not even light can escape from it, this is however a non-relativistic concept. The speed of 
light is not privileged in pre-relativistic physics, and a moving observer would not necessarily 
see it as a dark object. A more complete account of this history is given by Hawking and Ellis 
[39| including reprints of the original articles. 

The history of black holes proper dates back to just after the discovery of general relativity. 
The first non-trivial exact solution to the Einstein equations discovered by Schwarzschild 
in 1916 and named after him was, in fact, a black hole. It was however more than four 
decades before its properties were fully appreciated. The Kruskal-Szekeres extension of the 
Schwarzschild solution was discovered only in 1960. This was shortly followed by the discovery 
of the Kerr solution in 1963 representing spinning black holes. Its global properties were 
explained by Carter in 1966. The Kerr-Newman solution representing charged, spinning black 
holes was discovered in 1965. It was in 1964 that the phrase "black hole" was first coined 
by John Wheeler. During the same time, there were seminal developments in understanding 
the general properties of black holes beyond specific examples. This includes the study of the 
global properties of black hole spacetimes, the definition of event horizons, and crucially for 
the developments to be discussed here, the singularity theorems of Penrose and Hawking and 
the introduction of trapped surfaces by Penrose. This was soon followed by the understanding 
of black hole thermodynamics by Bekenstein, Bardeen, Carter and Hawking in 1973, and the 
discovery of Hawking radiation in 1974. The cosmic censorship hypothesis was formulated by 



1 



Penrose in 1979. The question of whether this is valid, i.e. if every singularity that results from 
the future evolution of generic, regular initial conditions is hidden behind an event horizon, 
is not settled and is one of the key unsolved questions in classical general relativity. The 
black hole uniqueness theorems which showed that the Kerr-Newman solutions are the unique 
globally stationary black hole solutions in Einstein-Maxwell theory in four dimensions was 
established in the 1980s following the work of Israel, Carter and Robinson. 

More recently, black holes have been the subject of intense study in quantum gravity where 
the calculation of black hole entropy has been seen as a key milestone for string theory and loop 
quantum gravity. There have also been important developments on the classical side where the 
long standing problem of calculating the gravitational wave signal from the the merger of two 
black holes was finally solved numerically in 2005 by Pretorius. In an astrophysical context, 
black holes are believed to be engines for some of the most violent events in out universe, such 
as active galactic nuclei. Astronomers have succeeded in locating a large number of black hole 
candidates with masses ranging from a few to billions of solar masses, and the direct detection 
of gravitational waves from binary black hole systems is expected later this decade. 

Most of these seminal developments have relied on event horizons to characterize the bound- 
ary of the black hole region (the singularity theorems are a notable exception). This is com- 
pletely reasonable when we are dealing with stationary situations, but can lead us astray in 
dynamical situations. As we shall elaborate later, event horizons are global notions and it is 
in principle not possible for a mortal observer to locate them. One possible alternative is to 
use the notion of trapped surfaces introduced by Penrose. While not entirely local since they 
are closed spacelike surfaces, these provide a quasi-local alternative which an observer could in 
principle locate in order to detect the presence of a black hole. Trapped surfaces lead logically 
to various kinds of quasi-local horizons including isolated, dynamical and trapping horizons. 
The goal of this chapter is to motivate and explain the quasi-local approach to studying black 
hole horizons, and to review some recent results. Somewhat surprisingly, we shall see there is 
still a major gap in our understanding of classical black holes in dynamical spacetimes. If we 
accept black holes as bonafide astrophysical objects, we still do not have a satisfactory notion 
of what the surface of a black hole is. Event horizons are not satisfactory because of their 
global properties, but there is as yet no established quasi-local alternative. 

The style of this chapter will typically be to start informally with simple examples and 
to use them as guidance for developing general concepts and definitions. In Sec. [2] we shall 
start with the simplest black hole, i.e. the spherically symmetric Schwarzschild spacetime, 
and understand the properties of its black hole region. This naturally motivates the funda- 
mental notions of event horizons and trapped surfaces, and to the boundary of the trapped 
region. As we shall see, the different reasonable definitions of the black hole horizon agree 
in Schwarzschild. This will not be the case in more general situations. Perhaps the simplest 
example is the imploding Vaidya spacetime which shall be our second example in Sec. |2.2| We 
shall see that at least in this simple spherically symmetric example, the location of the trapped 
surfaces can be determined. These two examples are then followed by general definitions of 
event horizons, and trapped surfaces in Sec. [S] which formalizes many notions introduced in 
Sec. [2] Some properties of trapped surfaces under deformations and time evolution are then 
discussed in Sec. |3.3| and this leads naturally to the notions of marginally trapped tubes, and 
trapping and dynamical horizons. We then restrict our attention to isolated horizons which 
describes the equilibrium case, when no matter or radiation is falling into the black hole (but 
the rest of spacetime is allowed to be dynamical). This is now well understood and we re- 
view the general formalism in Sec. [4] with the Kerr black hole as the prototypical example. 
In particular, we discuss two applications: black hole thermodynamics and the spacetime in 
the neighborhood of an isolated black hole. Sec. [S] reviews some results and applications for 
dynamical horizons and finally Sec. [6] provides a summary and some open issues. 
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The discussion of this chapter will be mostly self-contained, though digressions into the 
relevant mathematical concepts will be necessarily brief. Useful references for black holes 
and general relativity are [261 165] , and for more mathematically inclined readers [SD] is rec- 
ommended as a useful introduction to the relevant concepts in differential geometry. The 
discussion is meant to be accessible to physics graduate students who have taken a first course 
in general relativity (covering, say, the first part of the textbook by Wald [SS]). In the same 
spirit, the list of references is not meant to be exhaustive in any sense, and is mostly biased 
towards reviews and pedagogical material. The selection of topics is not meant to be exhaus- 
tive; this contribution is not to be viewed as a broad review article. It is rather a combination 
of pedagogically useful examples, and a brief description of some recent results. This mate- 
rial will hopefully whet the reader's appetite and motivate him/her to delve further into the 
subject. 

Some words on notation are in order. A spacetime is a smooth four-dimensional manifold 

M with a Lorentzian metric gab with signature ( h -|--|-). We shall use a combination of 

index-free notation and Penrose's abstract index notation for tensors lower-case Latin 
letters a,b,c, . . . will denote spacetime indices. Symmetrization of indices will be denoted by 
round brackets, e.g. X(ab) {Xat + Xba)/^, and anti-symmetrization by square brackets: 
X[ab] {Xab — Xba)/'2. The dcrivativc-operator compatible with Qab will be denoted Va, 
and the Riemnann tensor Rabcd will be defined by 2V[t,Vt,]Xc — Rabd^Xa for an arbitrary 
1-form Xa- The Ricci tensor and scalar are respectively Rab = Racb'^ and R = g°'^Rab- The 
coordinate derivative operator will be denoted d. The exterior derivative will be either denoted 
by indices, such as \/[aXt], or in index free notation as dX. The Lie derivative of an arbitrary 
tensor field T along a vector field X will be denoted ^xT. Where no confusion is likely to 
arise, we shall often not explicitly include the indices in geometric quantities. Unless otherwise 
mentioned, we shall work in geometrical units with G = 1 and c= \. We shall often deal with 
sub-manifolds of A^; a sub-manifold of unit co-dimension will be called a hyper-surface while 
lower dimensional manifolds (typically these will be 2-spheres topologically) will be called 
surfaces. All sub-manifolds shall be assumed to be sufficiently smooth. Unless mentioned 
otherwise, we shall be working with standard general relativity in four spacetime dimensions. 

2 Simple examples 

2.1 The trapped region in Schwarzschild spacetime 

We shall start by studying the gravitational field in the vicinity of a time-independent massive 
spherically-symmetric body. In this section we recall some basic properties of the Schwarzschild 
solution. 

The Schwarzschild metric is a static, spherically symmetric solution of the vacuum Einstein 
equations Rah = 0. It is usually presented as 

ds^ = - (^1 - ^ dt^ + (^1 _ dr" + v'idB' + sin^ d d<p^) . (1) 

Here r is a radial coordinate such that the area of spheres at fixed r and t is 47rr^; these spheres 
can be obtained invariantly by applying rotational isometries to a given initial point in the 
manifold. Each of these spheres are isometric to the standard round spheres in Euclidean 
space, and {9, (j>) are the usual polar coordinates. The time coordinate is t, and the metric 
is explicitly time independent in these coordinates. The parameter M is the mass and in 
non-geometrical units, we would have the replacement AI -> GM/(? . This metric turns out to 
be an excellent approximation to, say, the gravitational field in our solar system with the sun 
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treated as a point mass and with (f){r) — GM/(?r being its Newtonian gravitational potential. 
The quantity Rs ~ 2GM/c? is known as the Schwarzschild radius, and for the sun Rs ~ 3 km 
(which agrees with the Michell-Laplace idea mentioned at the very beginning of this chapter). 

The metric as written above is regular and non-degenerate for 2M < r < oo and — oo < 
t < oo. The singularity at r = 2M is not a true physical singularity [65] and can be removed 
by the transformation {t, r) — > («, r) where 

dv = dt+(l- — ]dr. (2) 



r J 

In these coordinates (the ingoing Eddington-Finkelstein coordinates) the metric becomes 

1 - du^ + Idv dr + r^dO.^ , (3) 

where dSl^ := d9^ + sin^ OdijP' . The metric is now regular for r > and — oo < ?; < oo. We 
could extend the solution further by going to double null coordinates, but this shall suffice for 
now. 

Consider now the vector fields 

— + -(^1-— ^— n--— U) 
dv 2 \ r J dr ' dr 

It is easy to check that these (future directed) vector fields are both null, i.e. £ ■ £ — n ■ n — 0, 
and £-n — —1. By convention, we take i to be outward pointing and n to be inward pointing; 
we have designated r — >■ oo to be "outwards" . 

Let us pause to recall the notions of expansion and shear of a vector field. For a timelike 
vector field the set of vectors orthogonal to ^ form a three-dimensional plane at each point. 
If A is an infinitesimal area element in this plane and A the affine parameter along ^, then the 
expansion of ^ is defined as 

Since a null vector field is orthogonal to itself, it is easy to show that any vector field V 
satisfying V ■ (, — can be written as a linear combination V — + /3e(i) -|- 7e(2) where eji) 
and 6(2) are mutually orthogonal unit spacelike vectors orthogonal to ^, and a, j3, 7 are real 
numbers. The expansion of ^ is then defined as in Eq. ^ above except that the relevant area 
element A is in the two-dimensional plane spanned by e(i) and 6(2). An alternate expression 
for the expansion which is usually more useful is 

9(0 = Q^'V^^ (6) 

where g"' is the (inverse of) the Riemannian metric in the (e(i),e(2)) plane. The trace of Va^t 
after projection is the expansion. The symmetric trace-free part and the antisymmetric parts 
give the shear aat and twist Uat respectively: 

(Jab = (V(aC6))^ - ^9(«)9a6 , UJab = (VfaCfc])^ ■ (7) 

where the symbol (• • • )^ indicates a projection in the (e(i), 6(2)) plane. The operators aab and 
LJab are responsible for transforming vectors in the (e(i),e(2)) plane. Consider now a set of 
neighboring null geodesies generated by Let be a connecting vector, i.e. it is transverse 
to r and is Lie dragged along = [C, C]" = ^'"^bC - C'^bC = 0. We get the effect 
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Figure 1: From left to right, an illustration of expansion, shear and twist respectively in the 
6(2)) plane transverse to as illustrated by the effect on circles in this plane. The expansion 
is an isotropic expansion, the shear is an expansion and contraction in orthogonal directions 

with the area being preserved, and the twist is a rotation. 



of the expansion, shear and twist as operators in the (e(i), 6(2)) plane leading to the evolution 
of in time 

The effect of expansion, shear and twist on is illustrated in Fig. [l] A particularly important 
result is the Raychaudhuri equation which gives the time derivative of the expansion for affinely 
parameterized nuU-geodesics: 

'^®(?) 1 r^2 ab , ab rr f a ^ 6 

~ ^ '^'^bCy +UJabi^ -Jab? t ■ (9) 



dX 



This is a particular component of the Einstein field equations and a derivation and applications 
can be found in e.g. |39l I65| . We shall have occasion to use it at various points during the 
course of this chapter. 

Returning to the vector fields i?" and defined in Eq. Q, we see that they are manifestly 
orthogonal to the constant {v,r) spheres. Thus their expansions must involve area elements 
on these spheres and it is in fact easy to calculate their expansions: 

r — 2M 2 

= e(„)(r-) = --. (10) 

For spheres outside the black hole region, i.e. spheres with r > 2M , we see that > 
and 0{n) < 0. This is how a round sphere in flat space behaves. However, for spheres in 
the black hole region, we get both expansions to be negative. Such spheres are known as 
trapped surfaces and play a fundamental role in black hole theory and, in particular, in the 
singularity theorems. The spheres on the r = 2M hyper-surface have 0(£) — 0, 0(n) < and 
are called margtnally trapped surfaces. Thus, we see that the r — 2M hyper-surface separates 
the region where the spherically symmetric trapped surfaces live and are a signature of a black 
hole spacetime. It is worth noting that the presence of trapped surfaces is not necessarily a 
signature of strong field gravity. The Riemann tensor (and thus the tidal force) is proportional 
to M/r^ . Thus at r = 2M, it is oc and large black holes have a correspondingly weaker 

curvature at their Schwarzschild radius. In this sense, trapped surfaces are a non-perturbative 
phenomena in general relativity. An observer falling into a sufficiently large black will not 
notice anything out of the ordinary. 
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The role of the r = 2M hyper-surface as the boundary of the region containing trapped 
surfaces (known as the trapped region) is typically not emphasized in standard textbooks on 
the subject. What is emphasized is instead the fact that starting from a point with r < 2M, 
there exist no timelike or null curves which can cross the r = 2M hyper-surface. This is easiest 
to visualize in a Penrose-Carter conformal diagram shown in Fig. [2] This is a convenient way of 
visualizing the r-t part of the Schwarzschild metric. Just as we had extended the Schwarzschild 
metric in Eq. ([T]) across r — 2M by using ingoing null coordinates, we can extend the metric 
of Eq. (|3| further by going to double null coordinates (u, v) where 

1 ] dr, dv = dt + ll 



du = dt-\^l jdr, dv = dt+ 1^1 — j '^^ ■ (^1) 

Fig. [2] is then obtained by performing a further re-scaling of coordinates and a conformal 
transformation which brings infinity to a finite distance. Details can be found in [65. 62] or 
in other standard textbooks on the subject. The original Schwarzschild metric is region I in 
this diagram while Eq. ([3| corresponds to I and II. Regions III and IV are mirror images of 
I and II respectively. In this figure, null curves are straight lines at 45°, is spatial infinity, 
i"*" is future timelike infinity and i~ is past timelike infinity. Future directed null curves in 
this figure end up either at the future singularity at r = (marked with a dashed line) or at 
future null-infinity labeled as J^""" (J^~ is past null-infinity). It is also worth pointing out that 
if region I is defined to be the "outside world" so that £°' is the outward pointing null normal 
(this is merely a matter of convention), then 0(f) < in regions II and III. However, region 

III has 0{„) > so that only region II has both expansions negative. Similarly, only region 

IV has both expansions positive. 

It is then clear that no point in region II can be connected with region I (or III) by a causal 
curve and thus justifies the term "black hole" for region II. The boundary of this region occurs 
at r = 2M and is called an event horizon. Note that the boundary of region IV also occurs 
at r = 2M but region IV is a time-reversed version of II; it has the property that no future 
directed causal curve can stay within it. Thus, region IV is called a white-hole. In the rest 
of this chapter, we shall only consider the portion of the r — 2M hyper-surface which bounds 
the black hole, i.e. region II. In a physical gravitational collapse situation, the white-hole (and 
region III) does not actually exist and is covered up by the matter fields which constitute the 
star; this will soon be seen explicitly when we study the Vaidya metric. 

We thus have two different routes for describing a black hole: trapped surfaces versus 
event horizons (though one might suspect a-priori that the two might be intimately related). 
Trapped surfaces seem to be more local than event horizons. To know that a particular null 
geodesic will not leave a particular region of spacetime, one might need to know properties of 
the spacetime far away from the starting point of the geodesic. On the other hand, for trapped 
surfaces, in the Schwarzschild case we have just needed the computations of the expansions 
in Eq. ( |10[ ) at a fixed value of {v,r). However, one shouldn't forget that the computations 
of the expansions are not at just a single spacetime point, but are instead to be performed 
at all points over a sphere. This issue is irrelevant for spherically symmetric trapped surfaces 
in spherically symmetric spacetimes, but it is an important point that one cannot identify a 
trapped surface by examining only a part of it. For this reason, the trapped surface condition 
is said to be quasi-local. In any case, for Schwarzschild, the two descriptions of the black hole 
region agree: the r — 2M hyper-surface is both the event horizon and also the boundary of 
the region of spacetime which contains trapped surfaces. 

Visualizing non-spherically symmetric trapped surfaces is harder, even in a spacetime as 
simple as Schwarzschild. As in many numerical investigations in general relativity, let us try 
to locate such surfaces on three-dimensional spatial hyper-surfaces. The equation 0(£) = 
turns into a minimization problem, and to a second-order elliptic equation in three dimensional 
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Figure 2: Penrose-Carter conformal diagram for the extended Schwarzschild spacetime. See text 
for details. 



space (we shall have more to say on this matter later) . The spatial hyper-surfaces depicted in 
Fig. |2] were all spherically symmetric, and thus a single curve in the Penrose diagram suffices 
for them. However, if we wish to depict non-spherically symmetric hyper-surfaces, we will 
need a collection of such curves, say one for each value of {6,4>)- An example is shown in 
Fig. [3] In this example, the spatial hyper-surface starts from on region III (but this detail 
is not important for our purposes and we could as well have started from on the left edge of 
region III). What is important is that the spatial hyper-surfaces, or alternatively all the curves 
shown in Fig. |3] intersect the event horizon. As we shall see later, due to the somewhat non- 
intuitive properties of such null surfaces, the intersection of such a spatial hyper-surface with 
the r = 2M hyper-surface is still a marginally trapped surface though now a non-symmetric 
one. More specifically, it turns out that £a is covariantly constant on the r — 2M surface so 
that Va^b, projected onto the r = 2M surface, vanishes identically (see the discussion around 
Eq. (561 and in Sec. 4.4 1. This means that all closed cross-sections of the r = 2M surface are 
marginally trapped. 

There would of course generally be non-spherically symmetric trapped surfaces on these 
spatial hyper-surfaces lying inside the marginally trapped one. Each spherically symmetric 
trapped and marginally trapped surface can be found by such a procedure. This construction 
clearly shows that there are many more non-symmetric trapped surfaces than symmetric ones; 
each spherically symmetric hyper-surface can be deformed in an infinite number of ways and 
still contain trapped and marginally trapped surfaces. 

We note finally that it is possible to come up with examples where part of the spatial 
hyper-surface extends arbitrarily close to the future singularity, but part of it is still outside 
the black hole region so that its intersection with the event horizon is not a complete sphere. 
There would then exist no marginally trapped surfaces on such a spatial hyper-surface [66] . 
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Figure 3: A non-spherically symmetric spatial hyper-surface in the Schwarzschild spacetime de- 
picted as a set of curves. 

2.2 The Vaidya spacetime 

As we have just seen, for a Schwarzschild black hole, all the natural definitions of the surface 
of a black hole agree. Thus, the r = 2M surface is both the boundary of the trapped region 
and also the event horizon. This has led to a widespread belief that the notion of a black hole 
and its surface is unambiguous. Matters are however not so simple in dynamical situations. 
Let us now look at what is perhaps the simplest example of a dynamical black hole, namely 
the spherically symmetric Vaidya spacetime [64] . 

The Vaidya solution is obtained by starting with the Schwarzschild metric in ingoing 
Eddington-Finkelstein coordinates of Eq. Q , and replacing the constant M by a non-decreasing 
function A'I(v): 




(12) 



The stress energy tensor for this metric is 



Tab ~ 



M(v) 



VaVVbV , 



where 



M{v) : 



dM{v) 
dv 



(13) 



This represents the collapse of null-dust to form a black hole and if M{v) > then Tab satisfies 
the dominant energy condition. If we choose a mass-function which is non-zero only for v > vq 
and constant after v = vi, then we will have (portions of) Minkowski and Schwarzschild 
spacetimes for v < vo and v > vi respectively. The Penrose-Carter diagram for this spacetime 
is shown in Fig. |4] A suitable set of null normals orthogonal to the constant (r, v) spheres are 





d 



(14) 



and their expansions are respectively found to be 




2 



r 



(15) 
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Thus, in this case, there are no spherically symmetric trapped surfaces outside the r = 2M[v) 
surfaces and as in Schwarzschild, the spheres with r — 2M{v) and fixed v are marginally 
trapped surfaces. 

The event horizon is also not difficult to locate. Consider the outgoing null geodesies 
generated by the vector field £ above. Some of these geodesies will reach infinity while others 
will terminate at the singularity. The event horizon is the boundary between the two cases. If 
we assume that the mass function reaches a finite final steady state value Moo, then the final 
black hole is a portion of Schwarzschild, and thus the condition r = 2Moo defines the final 
state of the event horizon. Thus, we want to find the outgoing null geodesic generated by i"' 



defined in Eq. (141 for which r — >■ 2Moo when v — >■ oo. Subject to this final state boundary 



rfr _ 1 / _ 2AfW . ^^^^ 



condition, we need to solve 

dv 2 

This is fairly easy to solve numerically for a generic mass function. 

Let us now note a few properties of this event horizon. A typical case is shown in Fig. [4] 
First note that the event horizon extends to the flat region. A mortal observer in the flat 
region who has no way of knowing that the gravitational collapse will occur at some time in 
the future, might actually be living near an event horizon. The existence of the event horizon 
has really no consequences for any physical experiment or observations that the observer can 
conduct locally, and contrary to popular belief, the observer can cross the event horizon without 
feeling anything out of the ordinary. Furthermore, even an observer in the intermediate region 
vo < V < vi, who can witness gravitational collapse occurring cannot know the true location 
of the event horizon. To illustrate this, consider Fig. [5] Here the mass function is non-zero for 
V < Vo as before, however, there are two phases. The mass function first reaches a constant 
value at vi but restarts again at a later time V2- The observer with v < V2 cannot know the 
value of Mao and can thus never know the true location of the event horizon. 

2.2.1 Examples of non-symmetric trapped surfaces 

In contrast to the event horizon, the r — 2M{v) surface seems to have the right properties. 
It bounds the region which has spherically symmetric trapped surfaces, it doesn't extend into 
the fiat region, it grows only when M > 0, it can be located quasi-locally, i.e. by checking the 
conditions for a trapped surface on a sphere, and it doesn't care about what happens to M{v) 
at late times. However, non-spherically symmetric trapped surfaces are not so well behaved. 
While there are no marginally trapped surfaces which lie completely within the flat region, we 
shall see that portions of them can extend into the flat region. 

One can try to find marginally trapped surfaces numerically. The standard procedure is to 
start with a particular spatial hyper-surface E, and to find a surface S in E for which 0(f) — 0. 
Let hab be the Riemannian metric on E induced by gab- If the unit normal on E is f", the 
unit spacelike normal to S on E is f° with = + f")/\/2 (see Fig. With this choice, 
noting that the metric on S is qab ~ hab — TaXb ~ gab + tatb — fafb, the condition = can 
be written as 

V2e(,) = ^q'^'Valb = Daf^ + Kabf'r' ~K = 0. (17) 
Here Kab = h'ahfW ctd is the extrinsic curvature, K is its trace, and D is the derivative operator 



on E (we shall explain these concepts is more detail later in Sec. 3.2 I. Taking coordinates 
{r,9,(j)) on S and assuming that the surface is given by an equation r = f{0,(j>), the above 
equation becomes a non-linear second-order partial differential equation for / which can be 
solved numerically. Typical methods assume that the surface is star-shaped, i.e. every ray 
from the origin r = intersects the surface exactly once; for a more complete description of 
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Figure 4: Penrose-Carter conformal diagram for the Vaidya spacetime. The region shaded in blue 
is flat and the region in red is isomorphic to a portion of Schwarzschild. The event horizon is 
labeled SH and is seen to be distinct from the r = 2M{v) surface. The two agree only in the final 
Schwarzschild portion. 
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Figure 5: Another Penrose-Carter conformal diagram for the Vaidya spacetime. The region shaded 
in blue is again the flat region. The mass function here has two phases where it is increasing. An 
observer with v < V2 could in principle locate the solid line marked 8%°^'^ but will not know that 
the mass function will increase again and that the true event horizon is in fact given by the surface 
(indicated by a dotted line in this figure) marked £%. 
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Figure 6: A closed marginally outer trapped surface on a spatial Cauchy surface E with intrinsic 
metric hab and extrinsic curvature Kab- The unit timelike normal to S is t"" and the outward unit 
spatial normal to S" is f°. A particular choice of the out- and in-going null normals are A={t°'±r^). 



this and other methods, we refer to [63]. We then choose a particular mass function and a E 
defined through a particular non-axisymmetric time coordinate and attempt to locate surfaces 
with 0(f) — 0. Let us review an illustrative result from a study reported in ^0] (see also [55] 
for another such study) The particular mass function chosen corresponds to a short pulse of 
radiation: 

, , for V <0, , , 

MovViv^ + W') for v>0. ^1^) 

The parameter Mo is the final mass and W determines the time-scale of the radiation pulse. 
We choose Mo — 1 and W = 0.1. The non-spherically symmetric time coordinate is taken to 
be 

t = V - r{l + acos9) . (19) 

The constant a determines the degree of asymmetry, and we choose q = 10/11. As in Fig. [s] 
spatial hyper-surfaces of constant t correspond to different sets of curves (one for each 0) in 
the Penrose-Carter diagram. Fig. [t] shows the marginal surface found on the t = —0.3 spatial 
hyper-surface. It shows the section in the x — z plane. The blue dotted line is the marginal 
surface, the green dashed line encloses the intersection of the fiat region with the hyper-surface, 
and the solid red curve shows the intersection with the r — 2M{v) surface, t = —0.3 surface 
(which is not a marginal surface). Thus, we see that the marginal surface extends in the 
fiat region, though in this example it is planar with 0(„) ~ (so it is not, strictly speaking, 
a marginally trapped surface). The marginal surface is seen to be only partially inside the 
r = 2M (v) surface. 



2.2.2 The trapped region 

Having looked at particular examples of trapped surfaces in a Vaidya spacetime, let us consider 
the trapped region, i.e. the portion of the manifold which contains trapped surfaces. The 
trapped region for Vaidya spacetime can in fact be studied analytically. The starting point 
for this goes back to a conjecture by Eardley in 1998 [30j : The outer boundary of the region 
containing outer trapped surfaces is the event horizon (an outer trapped surface has < 
and no restriction on 0(„)). One the one hand, it is known that trapped surfaces cannot cross 
the event horizon. On the other hand, in dynamical situations like Vaidya, the event horizon 
is growing in area and its cross-sections are not marginally trapped. Thus, while the outer 
trapped surface might get arbitrarily close to the event horizon, the limiting process is not 
trivial. The Vaidya spacetime provides a relatively simple setting to study this phenomena. 

Recent works in this direction have been |60U21|[^I23) V For any point p in the fiat portion 
of the black hole region of the Vaidya spacetime, Ben-Dov showed |21| that there exists an 
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Figure 7: A non-symmetric closed surface in a particular Vaidya spacetime with = and 
0(n) < 0. The solid red curve is the intersection of the spatial hyper-surface with the r = 2M{v) 
surface, the green dotted line is the boundary of the flat region, and the blue dashed line is the 
marginal surface located on this spatial hyper-surface. 



outer trapped surface S which contains p. This works even when p is arbitrarily close to the 
event horizon. In this case, most of the trapped surface actually lies inside the r = 2M{v) 
surface in the far future where v is large. There is a narrow "tendril" which is almost null for 
a large portion, and yet stretches from the far future right down to the fiat portion within the 
event horizon. This is precisely the kind of trapped surface whose existence was conjectured 
by Eardley in (30]. It is not clear that such highly non-symmetric trapped surfaces would be 
present in typical spatial hyper-surfaces used in numerical relativity simulations, or in fact, 
whether the standard numerical methods currently employed would be able to locate such a 
surface even if it were present. 

It is also possible to locate the boundary of closed future trapped surfaces (i.e. surfaces 
with Qf^i-) < and 0(„) < 0) in Vaidya spacetimes. The first result was obtained by Ben- 
Dov [5T], but the complete solution to the problem was found by Bengtsson & Senovilla 
[221 123) . Bengtsson & Senovilla have proved a number of results regarding the properties of 
trapped surfaces in spherically symmetric spacetimes, but here we shall only illustrate them 
by describing the past spacelike barrier for trapped surfaces in the Vaidya case. To this end, 
we need the following result (Theorem 4.1 of [23]): In a region 7?. of a spacetime, let 5" be a 
future pointing hyper-surface orthogonal vector field so that = —FVat for some F > and 
some i which increases to the future. If 5 is a future trapped surface which intersects TZ (but 
is not necessarily contained within TZ), then 5* cannot contain a local minimum of t at points 
with > 0. 

In a region TZ with a time coordinate such as t, the significance of this results is that once 
a future trapped surface enters such a region with initially decreasing i (i.e. if the surface is 
initially "bending downwards" in time), then t must continue decreasing. If TZ is bounded in 
the past by the event horizon, then clearly this result forces 5* to continue till it reaches the 
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event horizon. Since S cannot cross the event horizon (or even touch it), it becomes clear 
that the region TZ cannot contain even portions of future trapped surfaces which are bending 
downwards in time. In Vaidya, an appropriate ^ is the so-called Kodama vector 

Since ^a^" = —(1 — 2M{v)/r) it is clear that is future directed to the past of the r = 2M{v) 
surface. Furthermore, it is easy to check that 

^^gat = 2M(v)£alt ■ (21) 

Thus, q°'^J^^gab > if Af > 0. The surfaces of constant t, denoted by E* are spherically 
symmetric and defined by 

If M{v) is constant for v > v\, then there is a value of t = t such that E* coincides with the 
event horizon for v > v\. Alternatively, if M{v) asymptotes to a constant value, then there is 
a value oii = t such that E* also asymptotes to the event horizon. In both these cases, let 
us denote these E' as E. In the flat region, E* are horizontal lines in the Penrose diagrams. 
The behavior of E depends sensitively on M(v). In particular, it depends on the quantity 
:= lim„_»o M{v)/v. When p > 1/8, it is shown in [2^ th at S wi ll en ter the flat region, but 



in other cases it may not. These cases are depicted in Figs. 8a and 8b along with E for some 
other values of t. Any trapped surface which extends outside the r = 2M{v) surface must do 
so with increasing i and it must never "bend" downwards in time. If it does not do so, then 
the above result shows that it must continue downwards till it hits the event horizon where it 
must cease to be smooth if the condition 0(f) < condition is to be maintained (the event 
horizon is expanding and thus has positive expansion). It cannot terminate smoothly because 
this would imply the existence of a point on the trapped surface where t is a local minimum. 

A little thought then shows that E must then be a past barrier which future trapped 
surfaces cannot cross. Furthermore, it is also clear that there cannot be any compact future 
trapped surface contained entirely in the region to the past of E, and bounded between E 
and the event horizon; if there were, again there would have to be a point on the trapped 
surface where t is a local minimum. The boundary of the region containing trapped surfaces 
has thus been located. This boundary E is of course spherically symmetric (as one can prove 
on general grounds). However, it is not foliated by marginally trapped surfaces. Thus, as 
with the limit of outer trapped surfaces to the event horizon, it is clear that the the limit of 
marginally trapped surfaces to this boundary cannot be smooth. E is not a quasi-local object; 
it is as non-local as the horizon. Moreover, as far as we know, it does not have any features 
which might distinguish it as a black hole horizon. 



2.2.3 Lessons from spherical symmetry 

After this extensive discussion of the Schwarzschild and the Vaidya examples let us summarize 
the situation in spherical symmetry. We have seen the complications that can arise from 
having to consider non-spherically symmetric trapped surfaces in dynamical situations. In 
Schwarzschild there are no major surprises and any trapped surface can extend right up to the 
event horizon. This is not so in Vaidya; the obvious generalization of the r = 2M surface from 
Schwarzschild does not coincide with the event horizon. We need non-spherically symmetric 
trapped surfaces to "fill the gap" between the r = 2M{v) surface and the event horizon. It 
however turns out to be possible to study the trapped region in detail and to obtain a fairly 
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Figure 8: Surfaces of constant Kodama time for a Vaidya spacetime with /i := limy^o M (v) / v > 
1/8. In this case E lies partially in the flat region. A future trapped surface which extends in 
the flat region is depicted as S. It required to enter the r > 2M{v) region with increasing t. The 
right panel is similar to the left panel, but in this case E does not extend into the flat region. 
This happens for /i < 8. A similar picture also works when M(v) asymptotes to a finite value for 
V -> oo. This figure is essentially Fig. 15 of ,23j . 
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complete understanding of where trapped surfaces can (and cannot) occur. There turns out 
to be a difference between future- and outer-trapped surfaces (i.e. whether or not we consider 
the 0(n) < condition). Outer trapped surfaces can extend all the way to the event horizon, 
but understanding how such surfaces limit to the event horizon is subtle; there is a separate 
past barrier for future trapped surfaces which is distinct from the event horizon. 

What does this study tell us about non-spherically symmetric dynamical spacetimes? We 
note that the essential complication here, as noted by Eardley [30], is not that the black hole 
is non-spinning etc. Rather, the problem is to consider trapped surfaces which do not share 
the symmetry of the spacetime and to understand what happens to them as they are deformed 
towards the event horizon which will generally have positive expansion in dynamical situations. 
We might still expect outer trapped surfaces to extend to the event horizon in a similar fashion 
but there is, in general, probably no separate barrier for future trapped surfaces. 

The next step in this chapter will be to study how marginally trapped surfaces evolve in 
time and under general deformations, and this leads to the subject of quasi-local horizons. 
However, before doing so, we shall first formalize many of the ideas introduced in this section 
with some general definitions and results. 

3 General definitions and results: Trapped surfaces, 
stability and quasi-local horizons 

3.1 Event horizons 

The surface of a black hole is traditionally defined in terms of an event horizon which is 
the boundary of the region from where massive or mass-less particles can reach the outside 
world. The formal definition is however more involved, with the main difficulty being in how 
the "outside world" is to be defined. It is worthwhile to briefly sketch the various technical 
ingredients that go into the precise formal definition, if only to highlight once again the truly 
global nature of event horizons; see e.g. [651 139| for details. 

This requires one to attach future and past null-infinity J'^ as boundaries to the physical 
spacetime and to consider the causal past of J''^ . This causal past is the region of spacetime 
TL from which causal signals can escape to infinity and represents the "outside- world" . The 
future boundary of TL is the event horizon. In Fig. [2] the region TL is the union of regions I, 
II and III, the black hole is of course region II and the portion of the r — 2M surface which 
divides region II from I and III is the event horizon. A little thought shows that in order for 
this notion to capture the physical idea we have in mind, it is necessary to ensure that is 
complete in an appropriate sense. For example, if we were to look at the causal past of just 
a portion of J^^ even for Minkowski space, we would erroneously conclude the presence of 
a black hole in Minkowski space [33] • Similarly for Schwarzschild, we could end up with the 
wrong location of the event horizon if we looked at only a portion of J^'*^ (see again Fig. [2]). 

Since we need to construct a complete .^'^ and its global past in this definition, it is clear 
that the event horizon is a very global and teleological notion. There may in fact be an event 
horizon forming and growing in the room you are reading this article right now, because of 
possible events which might occur a billion years from now. An example is an observer in the 
flat region of the Vaidya spacetime in Fig. [4] This discussion also illustrates that formally, 
the notion of an event horizon cannot be used in a cosmological spacetime such as the one we 
inhabit, since it fails to be asymptotically fiat. 

In practice, the principle of calculating event horizons in numerical simulations is to start 
from an educated guess at late times, and to integrate a null geodesic or a null surface back- 
wards in time |63l [28] . The basis for these methods is the fact that when we integrate forwards 
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in time, a small initial error will cause the null geodesic or surface to diverge exponentially 
from the true solution and end up either in the singularity or at infinity. This implies that by 
integrating backwards from even a reasonable guess at late times one will converge to the true 
solution exponentially. 

3.2 Trapped surfaces 

Let's begin with the standard definitions of the I''' and 2"'^ fundamental form of a smooth 
non-degenerate sub-manifold A4 embedded in a spacetime {A4,gab)', our discussion mostly 
follows [50]. The fundamental form is just the induced metric metric hab on M so that for 
any vectors and Y'' tangent to M: 

habX^Y" ~ gabX^Y" . (23) 

If hab is non-degenerate, we can decompose the tangent space TpAf at any point p £ A4 as 

TpM = TpM e T^M , TpM n T^M = {0} , (24) 

where TpM is the tangent space to M and the subspace Tj^M is normal to it. Thus, an 
arbitrary non-vanishing vector field ^ defined at points of M can be split uniquely into a 
tangential and normal part: 

where • X = (25) 

for any vector X tangent to A4 . Then, for any X, Y tangent to A4 we have 

VxY = {VxYf + (VxY)^ . (26) 

The intrinsic covariant derivative on M is then defined as DxY ■- {VxY^. The 2"^^ funda- 
mental tensor 11 is an operator which takes two vectors tangent to A4 and produces a vector 
in the normal-subspace: TpM x TpM — ^ TpM: 

n(X,y) - (VxK)^. (27) 

It is east to show that D defined this way is a legitimate derivative operator, and that the 2"^"^ 
fundamental form is symmetric: n(X, Y) — Tl{Y, X). The symmetry of II is especially easy: 

n(x, Y) - n{Y, X) = (\/xY- VyX)^ = [x, y]^ = o . (28) 

In the last step we have used the Frobenius theorem which says that for a smooth sub-manifold 
M, if X, Y are tangent to M then so is their commutator. 

When M is a. hyper-surface, i.e. when it has co-dimension 1, then Np is 1-dimensional and 
spanned by the unit-normal A'^''. Thus we can define the extrinsic curvature Kab via 

n(X, Yy ■- -{KabX''Y'')N'' . (29) 

The symmetry of 11 implies that Kab = Kta- The most important case for us is however when 
the sub-manifold 5 is 2-dimensional and spacelike; the fundamental form, denoted by qab 
here, is a Riemannian metric on 5. We can again make the decomposition TpM = TpS (STjj^S . 
The normal space Tj^S is a l-f 1 dimensional Minkowski space. It will be convenient to choose 
two null vectors £ and n to span Tj^S. We are of course free to re-scale £ and n independently 
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by scalars, but we choose to use a normalization ^ ■ n = — 1 which cuts down the re-scahng 
freedom to 

e^Ai, n^A'^n. (30) 

We will always choose (£, n) to be future directed and i and n as outward and inward pointing 
respectively. The mean curvature vector K"' is the trace of the 2"^^ fundamental form: 

:= q'^'Hat' ^ -ein/' ~ OiDU . (31) 

The coefficients appearing here turn out to be precisely the expansions 0(e,n) discussed earlier. 



Under a re-scaling of the kind in Eq. (30 1, K" remains invariant. 

The different kinds of trapped surfaces correspond to properties of if". Two useful defini- 
tions are: 

• Future trapped surface {0(t) < 0, 0(n) < 0): -K"" is timelike and future-directed. 

• Marginally future trapped surface {0(t) = 0, 0{„) < 0): K'^ is null and future-directed. 

Furthermore, we shall consider only closed surfaces; this is an important condition as it, among 
other things, excludes trivial planar surfaces in flat space. 

It is also useful to remark on the physical significance of the condition 0(„) < 0. This 
condition holds for round spheres in fiat space and also, as we have seen, for round spheres 



in the trapped region of Schwarzschild (see Eq. (10 1. However, this may not be true for non- 
spherically symmetric trapped surfaces even in Schwarzschild. Furthermore, there are explicit 
black hole solutions found by Geroch and Hartle [33] which represent black holes distorted by 
surrounding matter fields which do not satisfy 0{n) < at all points on the horizon (but its 
average value over the horizon is still negative). In fact, for a number of key results that we 
shall mention below, it is not necessary to impose 0(„) < 0. Surfaces with vanishing outward 
expansion 0(f) — 0, with no restrictions on the sign of Q(„) will be called marginally outer 
trapped surfaces, usually abbreviated to MOTS. A surface with ©(^j < and no condition on 
0(„) will be said to be outer trapped. We have seen in the Vaidya example that there are 
differences in the location of trapped and outer trapped surfaces. 

The collection of closed future-trapped surfaces form the trapped region of the spacetime. 
More precisely, the trapped region 3^ consists of spacetime points which lie on a closed future- 
trapped surface. The boundary ^ oi 3^ is the trapping boundary. It is also common to 
consider these concepts restricted to a spacelike surface E, typically a Cauchy surface in an 
initial value problem set-up. The trapped region on E is the set of points on S which lie 
on a closed future-trapped surface contained entirely on E. The boundary of is denoted 

, and each connected component of is called an apparent horizon if it is the outermost 
such boundary on E. Since excludes trapped surfaces not contained in E, it is clear that 
SS^ G £^f]T.. It is important to not confuse the trapped region .fF' or its boundary .33 with the 



black hole region B defined in Sec. 3.1 and the event horizon. It can be shown [42j that with 
some additional regularity conditions, each connected component of the apparent horizon 
is actually a closed marginally trapped surface. The same result was proved in [5] with the 
regularity assumptions removed. In practice, this fact is what is used to locate the apparent 
horizon in numerical simulations. 



3.3 The stability of marginally trapped surfaces, trapping and 
dynamical horizons 

Let us now go beyond individual trapped/marginally trapped surfaces and look at their time 
evolutions. Consider a region of spacetime foliated by smooth spacelike surfaces E* depending 
on a real time parameter t. Start with initial data (the I''' and 2"^^ fundamental forms) at t — 
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Figure 9: The evolution of a MOTS in time. Start at t = with a MOTS 5*0 on a spatial hyper- 
surface 11°. Evolve the data on E° using the Einstein equations. If the evolution of the MOTS is 
smooth in time, then S will evolve to a MOTS St on S* at time t > 0. The collection of all the 
St will then form a smooth 3-surface H. We shall see that H will usually be spacelike so that the 
future null vectors and n° orthogonal to St (and the future null cone) will point "inwards" on 
T-l. The vector f" is the unit outwards normal to Sf on E*. 



and evolve it using the Einstein and matter field equations. This way, we obtain a solution 
to the field equations locally in time near E". The first question we wish to address is: If E° 
contains a MOTS 5*0, does it persist under time evolution and does it evolve smoothly? If 
it does evolve smoothly, then the union of all the MOTS 5* will form a smooth 3-surface H 
which we shall call a marginally trapped tube (MTT). A related question is then: How does H 
depend on the foliation E*? If we start with the same TP but choose the surfaces differently 
for t > (still requiring E* to form a smooth foliation), then will we still end up with a smooth 
MTT H'? If it exists, is it different from H? In numerical simulations, it is found that the 
apparent horizon can evolve discontinuously. Can this be understood analytically? 

It turns out that the answer to the above questions are intimately connected with the 
stability of So with respect to variations on E". To this end, we need to define the notion of 
the geometric vanation of a 2-surface 5* which is embedded in a spacetime M |54|[2ll4ll3]. Such 
a variation is a very general concept; it includes evolving in time, and an evolution following 
Einstein equations is a particular case. A smooth variation of a sub-manifold S is defined as 
a 1-parameter family of surfaces S\ (where A is a real parameter and takes values in some 
interval (— e, e)) such that: 5*0 is identical to S, each S\ is a smooth surface, and each point on 
S moves on a smooth curve as A is varied. We can then define a vector field as the tangent 
to these curves. This is depicted in Fig. [lO] We could also perform the variation along null 
directions or spatial directions in a given spacetime and of course, the variations do not need 
to be uniform on 5 and different points on S can move at very different speeds depending on 

If we have a relevant geometric quantity on S, we can compute it on each and differen- 
tiate it with respect to A, and the derivative is called the variation of that geometric quantity. 
If ff\ is such a geometric quantity (e.g. the expansion Oj^) on each S\), then we shall define 
Sqff :— d\ff\\\=o- It is also important to keep in mind that for a function /, in general the 
variation is not linear: Sfgff ^ f^qff- 
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For an MOTS defined on a spatial liyper-surface E, the relevant variation is along f°, the 
unit normal to 5*. The stability of S is meant to capture the idea that that if S is deformed 
outwards, it becomes untrapped, and an inward deformation leads to < 0. No condition 
on 0(„) is assumed. More precisely, the MOTS S is said to be stably outermost if there exists 
a function / > on S such that SfrO(e) > 0. S is strictly stably outermost if in addition 
SfrQ(e) 7^ somewhere on S. 

With this background, we can state the following result [31 If So is a smooth MOTS on 
E", and So is strictly stably outermost on E'^, then So evolves smoothly into smooth MOTSs 
St on E* at time t at least for sufficiently small (but non-zero) t. Furthermore, the union of 
the St forms a smooth 3-surface which we shall call "H. This holds at least as long as the St 
continue to remain strictly stable outermost. In addition, if Gai,^"^' > somewhere on S or 
if 1°' has non-vanishing shear somewhere on S, then H is spacelike. More generally, if the null 
energy condition holds then H is either spacelike or null. 

These results answer several of the questions raised at the beginning of this sub-section. If 
we were to start with the same Eo but choose different ones at later times, the MTT would 
still exist at least for some small time interval since the stability condition still holds at t = 0. 
However, H' and H would not necessarily coincide. The jumps that are observed in numerical 
simulations are because of the outermost condition: while an individual MTT can continue 
to evolve smoothly, a new MTT can appear further outwards. While not all MOTSs satisfy 
the stability condition, in all the cases that the author is aware of in numerical simulations, 
the MTTs continue to evolve smoothly. This suggests that the mathematical results could be 
strengthened. 

Having seen that there is a physically interesting class of MOTSs which evolve smoothly 
in time, it is reasonable to impose additional conditions on an MTT to capture the fact that 
they are black hole horizons. The first is the notion of a trapping horizon defined by Hayward 
|411I43| in 1994: A future- outer-trapping horizon (FOTH) is a smooth 3-dimensional manifold 
H foliated by compact 2-manifolds St such that 
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1. The surfaces S axe marginally trapped surfaces {0(e) = 0, 0(„) < 0), 

2. The directional derivative of along n" vanishes: ^„0^i) < 0. 

Historically, this was the first systematic definition of a quasi-local horizon and it played a key 
role in spurring further developments. At the time, it was not known whether a MOTS would 
evolve smoothly in time and if trapped surfaces would limit smoothly to the boundary of the 
trapping region (as we have seen in the Vaidya example, they in fact do not and this possibility 
was recognized in [41]). Using this definition, Hayward showed that H would be generically 
spacelike and is null if and only if the matter flux and the shear of £°' vanish identically (the 
result from [31 H] quoted previously is stronger) . Hayward also showed that it was possible 
to assign a surface gravity to the black hole, and to have versions of the laws of black hole 
mechanics applicable to H. The second condition seems similar to the stability condition 
defined above, but in fact the two may not necessarily agree. A FOTH requires implicitly 
that 1°' and have been extended smoothly in a neighborhood of 5* and does not refer to the 
variation of S. 

A dynamical horizon [141 115| is defined similarly, but it is designed explicitly for the case 
when the horizon is spacelike: A smooth 3-dimensional manifold T-i is said to be a dynamical 
horizon (DH) if it is foliated by compact 2-dimensional manifolds St such that 

1. The surfaces St are marginally trapped (0(f) = 0, 0(„) < 0). 

2. T-i is spacelike. 

The first condition is the same as for a FOTH, but the second condition specifies a-priori 
the causal nature of "H without any additional conditions on fields transverse to H. Even 
though the definitions of a FOTH and a DH are similar, neither implies the other. We shall 
return to dynamical and trapping horizons later in Sec. [S] but before that, we shall discuss 
the equilibrium case when the MTT is null in some detail. 

4 The equilibrium case: Isolated horizons 

Having discussed some general properties of trapped surfaces and their evolutions, we shall 
now consider a special case that is nevertheless of significant interest for various applications, 
namely when the black hole is in equilibrium with its surroundings and there is no matter 
or radiation falling into it. Since the work of the late 1990s and the subsequent years, our 
understanding of this special case has matured and and we can now consider it be well under- 
stood. All the well known globally stationary solutions in 4-dimensions, i.e the Kerr-Newman 
black holes, are included in this analysis. Also included are the dynamical cases when the 
black hole is itself in equilibrium but when the time dependent fields are relatively far away 
from the black hole. An example is a system of binary black holes when the separation of the 
black holes is much larger than either of the masses; to a very good approximation each black 
hole can be treated as being in equilibrium with its surroundings locally. In the same system, 
once the two black holes have coalesced, the final black hole will soon reach equilibrium after 
it's ringdown phase. The assumption of global stationarity obviously does not hold in these 
cases. Each of these cases is well modeled by the framework of isolated horizons that we shall 
now describe. We shall start with a prototypical example namely, the structure of the horizon 
of a Kerr black hole. This shall be followed by a general definitions and a summary of some 
key results. We shall conclude with two applications: the near horizon geometry and the 
mechanics of isolated horizons. It will be convenient to use the Newman-Penrose formalism 
in this discussion, and therefore we start with a short digression to discuss this formalism. 

We note that strictly speaking, it is not essential to use this formalism for many of the 
results that we will discuss later. However, it does prove to be very useful in a number of 
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cases. For our purposes, we discuss it here because it makes calculations very explicit and 
is thus useful for pedagogically; there are no tensor indices and all geometric quantities are 
expressed in terms of scalar functions which have clear geometric and physical interpretations. 
Since we wish to start with the Kerr black hole as an example, this formalism is particularly 
useful for studying its properties explicitly. 



4.1 The Newman-Penrose formalism 



The Newman-Penrose formalism [521 153] is a tetrad formalism where the tetrad elements are 
null vectors, which makes it especially well suited for studying null surfaces. See fSBi i2Bl 162] 
for pedagogical treatments (note that these references take the spacetime metric to have a 

signature of (H ) which is different from ours). Start with a null tetrad {£,n,m,fh) 

where £ and n are real null vectors, m is a complex null vector and m its complex conjugate. 
The tetrad is such that £ ■ n = —1, m ■ fh = 1, with all other inner products vanishing. The 
spacetime metric is thus given by 



gab = —tarn — Uah + mainb + mairib . 
Directional derivatives along the basis vectors are denoted as 



(32) 



(33) 



We shall see that the Newman-Penrose formalism employs almost all the Greek symbols, and 
therefore often leads to conflicts in notation. An example is the symbol A which is used for 
both the directional derivative along n" and for the isolated horizon itself; we will have soon 
have other such examples. Hopefully the notation should be clear from the context. 

The components of the connection are encoded in 12 complex scalars, the spin coefficients, 
defined via the directional derivatives of the tetrad vectors: 



Dl 


— (e + e)£ — Rm — Km , 


(34a) 


Dn 


— — (e + e)n + vrm + vrm , 


(34b) 


Dm 


— %£ — Kn + (e — e)m , 


(34c) 


M 


~ (7 + 7)^ ^ ^ '''"'^ ) 


(34d) 


An 


= — (7 + 7)w + T^m + vm , 


(34e) 


Am 


— i>£ — rn + (7 — 7)m , 


(34f) 


S£ 


— {a + I3)£ — pm — am , 


(34g) 


5n 


— —{a + P)n + jim + Am , 


(34h) 


5m 


= \£— an + (/3 — ct)m , 


(34i) 


Sm 


— fl£ — pn + {a — I3)m . 


(34j) 



While this kind of expansion may seem to some like a backwards step to the days before efficient 
tensor notation was developed, it is actually very convenient in cases where null vectors and 
null surfaces are involved. The use of complex functions is also efficient because it cuts down 
the number of quantities by half. 

Many of the spin coefficients have a clear geometric meaning, 
implies that £' 



0. 



is geodesic if and only if k = 0. 
Similarly, from Eq. 



(34a I 



(34e 



First note that Eq, 
Furthermore, it will be afiinely parameterized 
we see that n" is geodesic if and only if = and 
it is affinely parameterized if 7 -I- 7 = 0. Important for us in particular are the coefficients 
p and a which are related to the expansion, shear and twist defined earlier in Sec. 12.1| Let 
t"' be tangent to a null geodesic congruence, let it be affinely parameterized, and let C^" be a 
connecting vector for the congruence (it is transverse to £"" and satisfies [£, Q =0). As in Eq.js] 
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we need to look at Va^i, projected in the (m, m) plane. Note that in this plane, the metric is 
q"'' = 2m'-"m''', the antisymmetric area form is = 2im^°'m!'\ and the 2-dimensional space of 
symmetric trace-free 2"'^-rank tensors is spanned by m"'m'' and fh,°'fn}'. From the definitions 
above, it is easy to show that: 

0(1) = g^'Va^ = m^Sla + = -2Rep , (35) 

m}°'fri^^a^b = Imp, rn'rn'V a^b ~ —t^ ■ (36) 

Thus, we see that — 2Rep is the expansion. Imp is related to the twist, and a is the shear of 
t. Similarly, real and imaginary parts of fj, gives the expansion and twist of n", while A yields 
its shear. It is also easy to verify that 

[m, m]" — {jl — ^)t' + (p — p)n^ + (a — P)m°' + (/? — a)m" . (37) 

Thus, using the Frobenius theorem, we see that m and m can be integrated to yield a smooth 
surface if l'^ and n" are twist free. Furthermore, since the projection of Sm is determined by 
/3 — Q, it is clear that this determines the connection, and thus the curvature of this surface. 

Since the null tetrad is typically not a coordinate basis, the above definitions of the spin 
coefficients lead to non-trivial commutation relations: 

' (AD -DA)/ =(e + e)A/-h(7 + 7)D/-(r + 7r)5/-(r + 7r)5/, (38a) 

{5D~DS)f =(a + l3-n)Df + K,Af-(p + e-e)Sf-a5f^ (38b) 

' {SA-A5)f =-uDf + {T-a-^)Af + {ii-j + j)5f + '\Sf, (38c) 

i5S-SS)f =m^t^)Df + {p-p)Af + {a^p)5f^{a-l3)Sf. (38d) 

The Weyl tensor Cabcd breaks down into 5 complex scalars 

*o = Cabcdt^'ratrn' , *i = Cabcdt^rntT/ , $2 = Cabcdtrnrn'n'^ , (39a) 
*3 = Cabcdi"'n'frfn'^ , *4 = Cabcd'm'nrn'n'^ . (39b) 

Similarly, the Ricci tensor is decomposed into 4 real and 3 complex scalars "l?ij: 

^m^\Rabti\ On = ^i?,6(rn'' + m''m'), ^22 ^\Rabnn\ A=^, (40a) 
$01 = ^Rabtm'' , $02 = ^iiatm^m" , $12 = ^iJabm^n" , = . (40b) 

We are allowed to make transformations of the null-tetrad while preserving their inner prod- 
ucts, thereby leading to a representation of the proper Lorentz group. The allowed transfor- 
mations are parameterized by two real parameters {A, tp) and two complex numbers (a, b) (a 
total of six real parameters in all, as expected) 

(i) Boosts: ^ -> 4^ , n-^A'^n, m m 

(ii) Spin rotations in the m — fh plane: m — > e"^m , £ ^ £ , n ^ n 

(iii) Null rotations around £: £ ^ £ , m ^ m + a£ , n ^ n + am + afh + \a\^£ 

(iv) Null rotations around n: n ^ n , m ^ rn + b£ , £ ^ £ + hm + hfh -f |6pn 

The transformations of the spin coefficients and curvature components under these transfor- 
mations are not difficult to work out. Again, we refer to [561 1261 152] for a more complete 
discussion. 
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The relation between the spin coefficients and the curvature components lead to the so 
called Newman-Penrose field equations which are a set of 16 complex first order differential 
equations. The Bianchi identities, V[aRbc]de ~ 0, are written explicitly as 8 complex equations 
involving both the Weyl and Ricci tensor components, and 3 real equations involving only Ricci 
tensor components. See 1561 126( 1621 for the full set of field equations and Bianchi identities. 



4.2 The Kerr spacetime in the Newman-Penrose formahsm 

As the prototypical example for an isolated horizon, we now describe the structure of the 
Kerr black hole horizon. This will also illustrate the utility of the Newman-Penrose formalism 
when dealing with null surfaces. A detailed study of the various intricate properties of the Kerr 
spacetime can be found in [26]. Here we shall be brief and focus on the essential properties of 
the horizon. 

The Kerr metric with mass AI and spin a is usually presented in textbooks as (this is 
however not the form that Kerr originally derived it) 



l^^^)dt-+ P'dr^ ^J^Ml^dt d4> + p^de- + 5^ (41) 



where 



= r-^ + cos^ e, A = - 2Mr + = (r^ + a^)p^ + 2a^Mr sin^ 6 . (42) 

This metric has two Killing vectors: a timelike one = (flu)", and a spacelike rotational 
one if"' — (9^)". Based on the behavior of the metric at large distances, one can assign a 
mass Moo ~ M and angular momentum J^o = aM to the spacetime. There are multiple ways 
to justify this. Because of the existence of the two Killing vectors, the clearest definition is 
through the so-called Komar integrals [48] based on the two Killing vectors (see also |65j V 
Moreover, again based on the behavior of the gravitational field at infinity, one can assign two 
sets of higher multipole moments Mk and Jk |33l 1381 ITS] which turn out to be fully determined 
by M and a : + iJk = M{ia)'' , fc = 2, 3, . . . 

As in the original Schwarzschild metric, there are coordinate singularities when A = 0. 
This happens when 

r^r± = M± \JhP - a2 . (43) 
These can be removed by a coordinate transformation (t, r, Q, (j)) — >■ {v, r, 9, ip): 

dv — dt -\ — — dr , dip — dcp — -^dr . (44) 

This yields the metric in («, r, 6, tp) candidate 

'^^^'^\ ^ 2 , J o • 2., , AaMrsm^e 

— — dv 4- Idv dr — la sm 6dr dip dv dp 

p2 J p2 

2 ,.2 , E2sin2 6i 2 1.^. 
+ p dO -\ dip . (45) 

The horizon is the 3-dimensional surface r = r+ which we shall denote A. The intrinsic metric 
qab on A in (v, 9, ip) coordinates is obtained by setting r = r+ and dropping the dr terms in 
this metric. Rearranging terms we get 

2 ■ 2 a 

qab = " ""^f (V„i; - Q.-^dip) {Vbv - n-^dif) + plVaOVbO , (46) 
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where p\ := r\ + cos^ 9 and Q = a/2Mr+ — a/{r\ + a?). It is easy to verify that this 
metric has signature (0 + +) with the degenerate direction being 

Thus, the null normal to A acquires an angular velocity term in the presence of spin. 

The cross sections of this manifold, i.e. the surfaces of constant v, are spheres with a 
Riemannian metric qab- The area of such a sphere is time independent: oa = 47r(r-+ + a^). It 
is easy to verify that if we choose a different coordinate v' with the cross-sections still being 
complete spacelike spheres, the area of each cross section is still aA- 

A suitable choice of the ingoing and outgoing future directed null vectors are 

/r^ + a^\ a dad A d , . 

n^a^-i ^ — , ^ V„ = — + — + — ^- — ^ — . (48) 

\ J or ov + o(p lir'^ + a-') or 

On A agrees with the null direction given in Eq. ( |47[ ). The other null vector is clearly 
null because the metric does not have a dr^ term, and the scalar factor in n" is chosen to 
ensure I ■ n = —1. The covariant versions are: 

7"^ + / ^ . 2 ^„ ^ „ A „ ^ Aa sin^ 9 „ 

(49) 



A suitable choice for m" is 



a sin 6* „ (r'^ + a'^) sin&„ i 



\7aV + ^ \/a<p + -^pVa9 , 



a a sine 9 1 did . , 

V2p dv sin 9 dip ^pd9 

Here we have defined p := r + iacos9, so that p^ = \p\'^. It is unfortunate that the notation 
can be confusing. For example the p^ used in the Kerr metric is not to be confused with the 
spin coefficient p. This should hopefully not cause confusion because the spin coefficient p will 
vanish identically, and unless mentioned otherwise, will refer to + sin^ 9. 

We can now compute the spin coefficients at A and for the moment we shall restrict our 
attention to those spin coefficients which are intrinsic to A, i.e. do not require any derivatives 
transverse to A. These are: e, k, tt, a, j3, p, a, p, X. As is typical in tetrad formalisms, we do 
not need to compute any Christoffel symbols in order to compute any of the spin coefficients; 



the exterior derivative suffices. Using the definitions of Eq. (34 1 we can write the exterior 
derivatives of £a,ma,na in terms of the exterior products of the basis vectors. The spin 
coefficients are then combinations of contractions of the exterior derivatives with the basis 
vectors. As an example, the acceleration of ^" and its value at the Kerr horizon is 



e + e — £ n ■ 2Vr„t:i,i — — — — -, . (51) 

''^ 2Mr+ 2M{M + ^/KP - a^) ^ ' 

The other spin coefficients at the horizon turn out to be: 

T y/2arsui9 ,_„, 
>z = (T = p = \ = v = Q,Ti = a + l3 = -- ^ (52) 

Some of these can be understood on general grounds. First, since is tangent to a smooth 
surface, it must be hyper-surface orthogonal which means that we must have Imp — 0. Fur- 
thermore, if a null vector is hyper-surface orthogonal, it can be shown to be tangent to a 
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geodesic. This implies k = 0. The important condition on physical grounds is that the ex- 
pansion of 1°' vanishes: Rep = 0. Thus, as we saw for a Schwarzschild black hole, the cross 
sections of A are marginally outer trapped surfaces. Now let us turn to the Weyl tensor. It 
can be shown that the only non-zero component is 

*2 = -7 , (53) 

(r — jacos 6')'^ 

Can we now extract from these results general properties of a null surface which should behave 
like a black hole horizon in equilibrium? Can we assign physical quantities such as surface 
gravity, mass, angular momentum and higher multipole moments? We shall answer these 
questions in the next sub-section. Regarding angular momentum, we note that the spin 
coefHcient which vanishes for a = is tt. In fact, from the values of the spin coefficient, we 
can check that for any tangent to the horizon, X"Vaf'' ~ X"'LOat^ where 

LOa ■■— — (e -I- e)na + Trrria + Tffha ■ (54) 

We shall see that the angular part of uja yields the angular momentum of the horizon. 

4.3 A primer on null-hyper-surfaces 

The fundamental geometric objects in the theory of isolated horizons are null surfaces. Let us 
therefore start with with a discussion of the geometry of null surfaces in a Lorentzian manifold. 
The horizon of a Kerr black hole is a special kind of null surface, namely an expansion and 
shear free null surface. Such null surfaces are rather special from a geometrical point of view 
as we shall now explain. 



3.2 



Consider a smooth sub-manifold A4 oi a spacetime {A4,gab)- As earlier in Sec. 
define the l°*-fundamental tensor of A4, i.e. the induced metric hat as the restriction of gab 
to M: habX'^Y'' ■.=jjabX''Y^ for any two-arbitrary vector fields and tangent to M. 
The sub-manifold Ai is said to be null when hab is degenerate. In the non-degenerate case, 
the ambient covariant derivative operator V induces a natural derivative operator T> on M. 
Furthermore, D is the unique derivative operator compatible with hab, i.e. T>ahbc = 0. 

Can we repeat the steps described in Sec . (|3.2[ ) for defining the fundamental forms and 
intrinsic connection on a null surface? When M is & null hyper-surface, we shall call a null 
normal if it is along the degenerate direction of hab, so that habt"" = 0. We thus encounter a 



problem in the very first step, i.e. in the decomposition of Eq. (241. The null normal P is also 
tangent to M so that TpM n T^M 7^ {0}. The other route to defining D, namely finding 
the unique derivative operator compatible with hab doesn't work either because a degenerate 
metric does not uniquely determine a derivative operator. We can still define an intrinsic 
derivative operator D if we pick a particular subspace of TpM transverse to M. Following 
|19| ■ we first pick a spacelike subspace Sp of TpAi; for the Kerr horizon, a natural choice would 
be vectors tangent to the spherical cross-sections of fixed v. There will be two 1-dimensional 
sub-spaces of null vectors orthogonal to Sp. One of them is NpM, the null direction of hab 
and the other will be transverse to A4 which we shall call NpM, and we can decompose TpM 
as TpM © NpM © NpM. Associated to a particular null normal I"' we shall pick a vector n" 
in the transverse direction by requiring that I ■ n = We can then decompose ^ as 

i = e +e where ;l • (55) 

Here a and P are scalars, ^ is in <Sp at each p. We can then define the intrinsic derivative 
operator as before: VxY — {W xY)^ . However, T) would depend on our choice of Sp and 
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unlike in the non-degenerate case, there is in general no natural canonical choice. There is 
however one case when this is not an issue, namely when the 2"'^-fundamental form vanishes, 
V xY is always tangential to M and there is no need to decompose V xY ■ This happens when: 

/a^VftF" = -X^Y^Vaib = (56) 

for any X", tangent to M. Thus, la is covariantly constant on the null surface. 

An alternative way to state the same result (emphasized in e.g. [S]) is: if we start with a 
vector tangent to A and parallel transport it using the spacetime derivative operator V along 
a curve lying on A, then the vector remains tangent to A. Since parallel transport of y along 
is defined by X'^V aY^ = 0, this alternative criteria is also equivalent to laX^VbY" = 0. 



We shall see that Eq. ( 56 1 is satisfied for the various kinds of horizons that we shall now define. 



4.4 Non-expanding, weakly-isolated and isolated horizons 

Having understood imll-surfaces, we are now ready to define different kinds of isolated horizons 
with increasingly stronger conditions. We shall start with the minimum set of conditions, 
namely a marginally trapped tube which is null, and no condition on the ingoing expansion. 
A smooth 3-dimensional null surface A is said to be a non-expanding horizon if: 

• A has topology x R, and if : S'^ x R — >■ S''^ is the natural projection, then zu~^{x) 
for any x € are null curves on A. 

• The expansion :— q'^'^W atb of any null normal P of A vanishes. 

• The Einstein field equations hold at A, and the matter stress-energy tensor Tab is such 
that for any future directed null-normal —T^H!' is future causal. 

We shall consider only null tetrads adapted to A such that, at the horizon, coincides with 
a null-normal to A. We shall also consider a foliation of the horizon by spacelike spheres Sv 
with V a coordinate on the horizon which is also an affine parameter along I: Civ = 1; S shall 
denote a generic spherical cross-section of A. Null rotations about P correspond to changing 
the foliation. 

This deceptively simple definition of a non-expanding horizon leads to a number of impor- 
tant results which we state here without proof, most of which are however well illustrated by 
the Kerr example discussed earlier. 

1. Any null normal P is a symmetry of the intrinsic degenerate metric qab on A: J^iqab ~ 0. 

2. The null normal of A is only given to be expansion free. However, a non-expanding 
horizon is also shear free. To show this, we use the Raychaudhuri equation: 

^eQi = KeQ^i^ - ^6^,, - - Rabti" . (57) 

Setting 0{^) — 0, and observing that the energy condition implies Rabi"'^'' > 0, we get 
that the sum of two non-negative quantities must vanish. 

\a\' + Rabre'' = 0. (58) 

This can only happen if a = and Rabt-°'(!' = on the horizon. Thus, the full projection 
of Va^t on A vanishes, and as we saw in Sec. |4.3| this is just the condition required to 
ensure that the induced derivative operator on A is well defined. 
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3. The Weyl tensor components and vanish on the horizon. This imphes that »I'2 is 
an invariant on A as long as the null-tetrad is adapted to the horizon; it is automatically 
invariant under boosts and spin rotations (it has spin weight 0), and it is invariant 
under null rotations around I because ^'o and vanish. Similarly, the Maxwell field 
component vanishes on the horizon, and is invariant on A. Both ^2 and (^\ are 
also time independent on the horizon. 

4. There exists a 1-form oja such that, for any vector field X°' tangent to A, 

X^Va^* = X^LJat ■ (59) 

The 1-form LOa plays a fundamental role in what follows. The puUback of lOa to the 
cross-sections S will be denoted cja- 

5. The surface gravity of I is 

K(^) = tuja ■ (60) 

We will say that A is extremal if K(^) = and non-extremal otherwise. Here we shall 
always assume that A is non-extremal. The curl and divergence of lo carry important 
physical information. The curl is related to the imaginary part of the Weyl tensor on 
the horizon 

do; = Im[*2]^e. (61) 
and its divergence specifies the foliation of A by spheres [^. 

6. By the geometry of A, we shall mean the pair (gab, Ca). Clearly, qab yields a Riemannian 
metric qab on the cross-sections of A. In turn, "Da is determined by lOa and by the unique 
derivative operator T>a compatible with qab- 

We need to strengthen the conditions of a non-expanding horizon for various physical situa- 
tions. The minimum extra condition required for black hole thermodynamics and to have a 
well defined action principle with A as an inner boundary of a portion of spacetime, is formu- 
lated as a weakly isolated horizon [12]: A weakly isolated horizon (A, is a non-expanding 
horizon equipped with an equivalence class of null normals \(\ related by constant positive 
rescalings and such that 

Li^^a = . (62) 

If we recsale I — > to a transforms as — > oj^ -f 9t, In /. It is thus invariant under constant 
rescalings and there is a unique lOa corresponding to the equivalence class [£]. 

The zeroth law holds on weakly isolated horizons, i.e. = CiOa is constant on A: 

^lUJa = ^"IV^b^a] + VaituJb) = ©aS(,) . (63) 



In the second step we have used Eq. (61 1 to conclude that i x 22?[5a;„] = Im^ [^'2] eta = 0. 
Note that under a re-scaling ^° fi , u}a transforms as cja — >■ tUa 4- Da In / so that it is 
invariant under constant rescalings. 

Any non-expanding horizon can be made into a weakly isolated horizon by suitably scaling 
the null generators. Thus, the restriction to weakly isolated horizons is not a genuine physical 
restriction. One could go ahead and impose further physical restrictions on the intrinsic horizon 
geometry by requiring that not only lOa, but the full connection Da on A is preserved by I"': 
An isolated horizon (A, [P\) is thus a non-expanding horizon equipped with a equivalence class 
of null- normals related by constant positive rescalings such that [Ci,T>\ = [9]. 

It can be shown that the gauge-invariant geometry of an axisymmetric isolated horizon can 
be fully specified by the area oa and two set of multipole moments M„, J„ for n = 0, 1, 2, . . . 
[10| . In contrast to the field multipole moments defined at infinity, M„ and J„ are the source 
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multipole moments. We shall not discuss these moments here in any detail, except to say 
that in the vacuum case, the moments are essentially obtained by decomposing at A into 
spherical harmonics based on preferred coordinates adapted to the axial symmetry. As at 
infinity, the Kerr horizon corresponds to a specific choice of these multipoles, but due to the 
non-linearity of general relativity, the field and source moments will not generally agree. See 
e.g. [611 1471 H51 159) for some applications of these multipole moments. 

4.5 The near horizon geometry 

In a number of astrophysical applications where black holes play a role, it is not directly 
the horizon which is involved but rather the spacetime in the vicinity of the black hole. An 
especially interesting example of relevance to gravitational wave observations is the case of a 
binary system consisting of a stellar mass black hole or star orbiting around a much larger 
black hole (see e.g. :1, for a review of the astrophysics of such systems). As the small particle 
orbits the large black hole, it effectively "maps" the spacetime in the vicinity of the large black 
hole, much as the motion of satellites around Earth enables us to map Earth's gravitational 
potential and therefore its shape. If carried out to sufficient precision, a measurement of 
gravitational waves from such a system would enable us to determine the gravitational field 
in the vicinity of the large black hole (see e.g. [581 144[ ) . One expects that the black hole is 
well approximated by a Kerr spacetime and we can seek to measure deviations from it thereby 
testing an important prediction of general relativity. Typical studies on this subject assume 
the black hole to be Kerr, which is entirely reasonable to a very good approximation (nearby 
stars or other matter fields might distort the black hole somewhat, but this is expected to be a 
small effect). However, for mathematical purposes, we could pose the question from a different 
viewpoint: If we specify the intrinsic geometry of the horizon, then to what extent is the near 
horizon spacetime determined by the horizon geometry? If we assume the large black hole to 
be in equilibrium (an excellent approximation for the case we have just described), then it 
seems reasonable to model the black hole as an isolated horizon. So the question then is: Can 
we find solutions to the Einstein field equations which admit a generic isolated horizon as an 
inner boundary? If so, then what is the extra data (beyond the intrinsic horizon geometry) 
that needs to be specified? (A full solution to the problem would require us to go beyond 
isolated horizons an to consider small deviations from equilibrium, but we shall not discuss 
this generalization here). 

It turns out that these questions can be clearly answered if we use the characteristic initial 
value formulation of Einstein's equations where free data is specified on a set of intersecting 
null hyper-surfaces [311 1571 1321 162| . Consider A'^ dependent variables "(/"/(I = 1, . . . , A'^) on a 
spacetime manifold with coordinates a;". We shall be concerned with hyperbolic first-order 
quasi-linear equations of the form 



In the standard Cauchy problem, one specifies the tpi at some initial time. A solution is then 
guaranteed to be unique and to exist at least locally in time. The characteristic formulation 
considers a pair of null surfaces Ao and A/i whose intersection is a co-dimension-2 spacelike 
surface S. It turns out to be possible to specify appropriate data on the null surfaces and on 
S such that the above system of equations is well posed and has a unique solution, at least 
locally near S. 

In our case, the appropriate free data is specified on the horizon and on an outgoing past 
light cone originating from a cross section of the horizon. Such a construction in the context of 




(64) 
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Figure 11: The near horizon coordinates. The isolated horizon is A and the transverse null surface 
is A^. The afhne parameter along the outgoing null geodesies on Af is r, and w is a coordinate along 
the null generators on A, and are coordinates on the cross-sections of A. 



isolated horizons was first studied by Lewandowski [51' who characterized the general solution 
of Einstein equations admitting an isolated horizon. This was worked out in detail in [49] 
which we follow here; similar results from a somewhat different perspective are discussed in 



|25| . The general scenario is sketched in Figure 11 We consider a portion of the horizon A 
which is isolated, in the sense that no matter and/or radiation is falling into this portion of 
the horizon. For a cross-section 5", the past-outgoing light cone is denoted by A/". The null 
generators of A and A/" are parameterized by v and r respectively; are coordinates on S. 
This leads to a coordinate system (u,r, a;') which is valid till the null geodesies on A/" start to 
cross. The field equations are solved in a power series in r away from the horizon. 

Let us now assume that the vacuum Einstein equations hold in a neighborhood of the 
horizon A. Following [11] we introduce a coordinate system and null tetrad in the vicinity of 
A analogous to the Bondi coordinates near null infinity. See Fig. jllj Choose a particular null 
normal i"' on A. Let v be the afiine parameter along so that £°''VaV = 1. Let 5*1, denote 
the spheres of constant v. Introduce coordinates a;' {i — 2, 3) on any one Si, (call this sphere 
So) and require them to be constant along i'^: £°'\/aX^ ~ 0; this leads to a coordinate system 
on A. Let be a future directed inward pointing null vector orthogonal to the Sv 
and normalized such that i ■ n = —1. Extend n" off A geodesically, with r being an afhne 
parameter along — n"; set r = at A. This yields a family of null surfaces A/'„ parameterized 
by V and orthogonal to the spheres S„. Set (VjX"^) to be constant along the integral curves 
of n" to obtain a coordinate system {v, r, x^) in a neighborhood of A. Choose a complex null 
vector tangent to Sq. Lie drag m" along l"": 

£.em'' = on A. (65) 

We thus obtain a null tetrad {£, n, m, m) on A. Finally, parallel transport I and m along — 
to obtain a null tetrad in the neighborhood of A. This construction is fixed up to the choice 
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of the and m" on an initial cross-section So- We are allowed to perform an arbitrary spin 



transformation m 



m on So- 



With the Bondi-like coordinate system in hand, we can now in principle use the coordinate 
basis vectors in the {v, r, a;') coordinates to construct an arbitrary null-tetrad near the horizon. 
The evolution equations for the component functions of the null tetrad will follow from the 

above construction. Let us start with ria and n". Wc have the family of null surfaces J^v 
parameterized by v; Ua is normal to the Afv, and r is an afHne parameter along — This 
implies that we can choose 



ria = -da 



and 



n"Va 



dr 



(66) 



To satisfy the inner-product relations i°'na = — 1 and m°'na = 0, the other basis vectors must 
be of the form: 



^ 9 d d ^ ^ d d 

■ a-D = — + U^+X^--, m'^Va--=6 = n-+e 



dv dr 



dr dx^ 



(67) 



The frame functions Z7, X'' are real while fi, (J' are complex. We wish to now specialize to the 
case when is a null normal of A so that the null tetrad is adapted to the horizon. Since 

dv is tangent to the null generators of A, this clearly requires that U,X^ must vanish on the 
horizon. Similarly, we want to be tangent to the spheres Sv at the horizon, so Q. should 
also vanish on A. Thus, f7, X',f2 are all 0{r) functions. 

To expand the metric in powers of r, we start with the frame fields, and the radial frame 
equations derived from the commutation relations. The strategy is the same as for the spin 
coefficients. The radial equations give us the first radial derivatives by substituting the horizon 
values on the right hand side, and taking higher derivatives leads to the higher order terms. 
The calculations are straightforward and lead to the following expansions: 



U = rk + r'^ { 2 7r'°M 



+ Re[4'^°^] ) + 0(r 



1 = 



X' = nc^o) + f2C(o) + ©('■'^ 



,(o)\(o) 



+ 2*4 



(0) 



4) + 0{t^ 



(68a) 

(68b) 
(68c) 

(68d) 



The contravariant metric is seen to be given in terms of the frame fields as follows: 



ff'-'- = 2([/ + |n|^), 5'"- = i, 



(69a) 
(69b) 



The null co-tetrad can be calculated easily up to 0{r^) 
n = —dv , 



(70) 



dr - («r + Re[*(°V') dv - (^tt^^V + r^) ^f^dx' - (^n^°^r + ^*i°V^) ^i°^dx' , (71) 



(0) 2 



r^]dv + {l- ii^°\)ifUx 



(A(°V+i* 



(72) 
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Here ^'^^ are defined by the relations Ci°''C(o) ~ ^^'^ Ci°'C(o) = 1; it will be convenient to set 



,(0) ._ AO) 



~ £.i daX^. In powers of r, the metric is: 

gat = -2£(„n,) + 2m(„m,) = g^°^ + g^,\ + ^gf.V^ + • • • , (73) 



where 



(75) 



5l",^=29,„r9,)t; + 2mJ"„^<;\ (74) 
+ 4M(°'m[»'.nl«) + 2A(»)m[>l») + 2A(»)^»',nl« 

+ 4 + A(°V(°) - *(°)) mf>,)« + 4 (^(/.f"')^ + |a(°)Q m[0'm(°) (76) 



Iterating this procedure to higher orders is, in principle, straightforward. This calculation 
provides a starting point for a number of ongoing work in applying isolated and quasi-local 
horizons to astrophysical situations. 



4.6 Angular momentum, mass, and the first law for isolated 
horizons 

The first law for black holes, and black hole thermodynamics in general, was developed by 
Bekenstein [20], and by Bardeen, Carter and Hawking |17| in 1973 in analogy with the laws of 
thermodynamics. The he zeroth law says that the surface gravity is constant over the black 
hole horizon. We have already seen that this is true for a weakly isolated horizon: the surface 
gravity Ke — £°'LJa is constant on A. The main difference with the standard formulation for a 
globally stationary the surface gravity refers to the globally defined stationary Killing vector 
which is normalized to have unit norm at infinity. A weakly isolated horizon does not refer to 
any globally defined Killing vector. The standard formulation of the second law says that the 
area of the event horizon can never decrease in time. The area of a non-expanding horizon is 
constant, thus the second law is trivial in this context. 

Let us now turn to the first law. The standard formulation for a stationary black hole is 

SM^^Sa + nSJ. (77) 

Here M is the mass measured at spatial infinity, k is the surface gravity at the event horizon 
but using a vector field normalized at infinity, a is the area of the horizon, Q is the angular 
velocity at the horizon, J is the angular momentum at infinity. Electromagnetic fields can 
also be included and leads to additional terms This is reasonable when applied to a globally 
stationary spacetime, but clearly needs to be refined for a black hole in equilibrium locally 
in an otherwise dynamical spacetime. It turns out that it is possible to formulate the first 
law for isolated horizons using quantities defined only at the horizon, without reference to the 
behavior of fields at infinity [G] [T] 1121 [8] . The set-up is a variational problem in a portion of 
spacetime bounded inside by an axisymmetric weakly isolated horizon (A, [i""]) (the extension 
to multiple horizons is straightforward). 

Before talking about the first law, we need to first have suitable notions of the horizon 
angular momentum and mass. We begin with angular momentum. Fix an axial symmetry 



32 



at the horizon. This means that tp" must preserve (i) the equivalence class [i""] of null normals 
that is prescribed for the weakly isolated horizon, (ii) the intrinsic metric qab and (iii) the 1- 
form LJa- Furthermore, (/p" should commute with and it should look like a rotational vector 
field in that it should have closed integral curves, an affine length of 2n, and should vanish at 
exactly two null horizon generators. Consider then a rotational vector field 0" in spacetime 
such that at A it is equal to the fixed symmetry: </>°|^ — f"". At infinity, we require that 
it approach some fixed rotational symmetry of the asymptotic flat metric. We then need to 
find the Hamiltoniaij^ which generates motions along 0". The Hamiltonian can be shown 
to reduce to integrals over the boundaries at the horizon and at infinity. The term at A is 
identified with the angular momentum of A: 

jf = (f ^°'^a'^<^- (78) 

Similarly, the notion of energy corresponds to time translations. Thus, we consider time 
evolution vector fields t"" on spacetime such that at the horizon, it approaches a general 
symmetry vector Ai"^ + Here the coefficients A and Q are constants on A, but are 

allowed to vary in phase space. The strategy is then again to compute the surface term at A 
in the Hamiltonian Ht which generates motions along f", and the surface term at A is to be 
identified as the energy Ef^. Surprisingly, it turns out that motions along t" are not always 
Hamiltonian and in fact, the Hamiltonian exists if and only if 

5E^ = ^SaA + ntSJA ■ (79) 

This is just the first law, and a vector field is said to be permissible if the first law for Ef^ 
is satisfied. If the first law is satisfied, it is easy to see that Kt,^t can depend only on the 
horizon quantities (oa, J a), and must satisfy integrability condition 

dnt{aA, Ja) _ d^{aA, Ja) 



OJa daA 



(80) 

This ensures that the right-hand-side of Eq. ( |79| l can be integrated to yield an exact variation, 
and to thus have a well defined Ef^ . 

A preferred choice of at the horizon can be obtained by choosing K{aA, Ja) and Q{aA, Ja) 
to have the same functional dependence on {a a, J a) as in the Kerr spacetime: 

Ra - 4Ja ^ 2 Ja 



This leads to the horizon mass: 



" = „ ..0 ■ (81) 



A^A = ^^i?i+4Ji. (82) 

Finally, we note that is independent of the choice of cross-section S and even though it 
requires a weakly isolated horizon to carry out the Hamiltonian computation, the formula of 



Eq. ( 78 1 itself is well defined even on a non-expanding horizon. If a cross-section 5 of A is 
contained within a spatial hyper-surface E, and if is the unit spacelike normal to S in E 
and Kab is the extrinsic curvature of E, then we can rewrite as 

Jl = ^ f> K^iV^^^- (83) 



^Recall that in a phase space, a Hamiltonian is responsible for generating time evolution via the Poisson bracket. 
Thus, for any function F in a phase space, F = {H, F}. In the present case, the phase space consists of gravitational 
(and other) fields which satisfy the appropriate boundary conditions. 
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This is particularly convenient in numerical relativity where one routinely located MTSs on 
spatial Cauchy surfaces, and would like to use them to characterize the properties of a black 
hole in real time while the simulation is in progress |29] . The computation of angular momen- 
tum is a surface integral over the MTS and the mass is just an algebraic expression. These 
methods are now in common use in numerical simulations. 



5 Dynamical horizons 

5.1 The area increase law 

The second law for event horizons states that the area of an event horizon can never decrease in 
time. This was first suggested by Bekenstein and is the starting point for associating the area 
of a black hole horizon with entropy. We have thus far not talked about the area increase law 
for quasi-local horizons except in the context of isolated horizons where it is essentially trivial. 
For a dynamical horizon, the area increase law follows easily from the fact that 0(n) < 0. Let 
"H be a dynamical horizon and S a generic MTS on it. Let be the spacelike unit normal to 
S within H; this is not to be confused with the unit normal f"' to S which lies on a spatial 
hyper-surface intersecting H as in Fig. [9] Let t'^ be the unit-timelike vector normal to H. We 
assume also that points outwards, in the sense that if f° is the outward normal to S on 
a spatial hyper-surface E, then r^r^ > 0. Finally, let H be bounded by the cross-sections 5*1 
and 52. Then, a suitable choice for the out- and in-going null normals to 5* are 

Then, if Qat is the intrinsic Riemannian metric on S: 

V2q'"'Varb = q^'Vait - r") = Gf^) - e(„) > . (85) 

Thus, the area element on S and thus its area increases along r". (Though we shall not pursue 
this further, one could imagine relaxing the requirement 6(n) < by an average on S and still 
obtain the same result). This is the area increase law for dynamical horizons. 

We can go further and ask whether it is possible to obtain a physical process version of 
the area increase law which relates the increase in area from an initial cross-section 5*1 to 
a later time 5*2 to the amount of energy or radiation falling into the black hole between 5*1 
and 82- An early result along these lines was proved by Hartle and Hawking in 1972 [40^ for 
perturbations of the event horizon. We note however that such a law does not exist for event 
horizons in general. A case in point being the Vaidya solution which as we have seen, grows 
in flat space when nothing falls into the black hole. Let us then consider the area increase law 
on a dynamical horizon H. Let us first consider angular momentum and the change in angular 
momentum from 5*1 to 52. 

Since "H is a spacelike hyper-surface we have, as discussed earlier, the induced metric hab, 
the associated derivative operator Da, and the extrinsic curvature Kab- As before, let 
be the unit timelike normal to H and let be the unit spacelike normal to a cross-section 
5 within H. The Einstein equations then show that (hat, Kab) cannot be specified freely, 
but must satisfy the Hamiltonian and momentum constraint equations (see e.g. [65])- The 
momentum constraint is: 

DbiK'''' - Kh"'') = 8nT'"'n,ht . (86) 
Contracting both sides with a rotational vector field and integrate by parts to obtain 

^ i Kab^^r'd^V -^i Kab^^r' = / fT,t,r>'' + -^P'^'^^hab) , (87) 
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where P'^'' = — Kh"'''. We then identify the angular momentum of a cross-section S as 



4"' / KabV^'^r' d'V . (88) 



s 



Eq. ( 87 1 is thus a balance law for angular momentum [141 [15] . It relates Jg^ — Jg^ with the 
gravitational and matter flux crossing the horizon between 5*2 and 5*1. We note that if i^" is 
a Killing vector of the metric hab, then the gravitational contribution vanishes identically. 

A similar (but more involved calculation) leads to a balance law for the change in area, or 
rather for the area radius R defined as i? = \/a/47r [141 115| : 



R2 

2 



= / r.f.r^C'rfV + / Nr(\a\^ + 2\(:f) d^V. (89) 



.|2 ^ ab 



Here \a\ := cTabO-" with aab being the shear of the outgoing null vector = r"" + r"; \(^\ 



(aC where C = h^^r^Vdb; and iV^ — \DaRD''R\^l'^ . Eq. |89| is the desired area balance 
law. Again the flux terms consist of a matter and gravitational contributions. The integrand 
in the gravitational part is manifestly local and non-negative, and the matter contribution is 
positive if the dominant energy condition holds. The gravitational contribution in fact vanishes 
in spherical symmetry as it should. Along similar lines, Jaramillo and Gourgoulhon obtained 
a second order differential equation for the area [37] which leads to a causal evolution of the 
area subject to initial conditions (the analogous evolution of the area of an event horizon [27] 
turns out to be teleological in the sense that it requires one to specify a boundary condition 
near future timelike infinity). 



Let us conclude with few words about angular momentum. We have identified Eq. (88 1 



with the angular momentum of a cross-section of the dynamical horizon, while earlier we 



had identified Eq. ( 83 \ with the angular momentum of a cross-section of an isolated horizon. 
It can be shown that when is a symmetry of the intrinsic metric, then the two agree. 
Recently, it was shown by Jaramillo et al |46) . that for an axisymmetric MOTS which is stably 
outermost, and if the spacetime satisfies the dominant energy condition, then the angular 
momentum J defined as above, satisfies the inequality |J| < a/Svr where a is the area of 
S. Here, axisymmetry is imposed only at S and not globally. This surprising result further 



validates the identification of Eq. ( 88 1 as the angular momentum. 



5.2 Uniqueness results for dynamical horizons 

In earlier sections, we have already discussed the location and time evolution of trapped 
and marginally trapped surfaces. We now discuss further results obtained by Ashtekar and 
Galloway [13] related to the issue of uniqueness of marginally trapped surfaces in the dynamical 
case. 

The first result, proved in 13), concerns the uniqueness of the foliation of a dynamical 
horizon by marginally trapped surfaces. Let "H be a dynamical horizon foliated by a set 
marginally trapped surfaces St with t being a continuous real parameter taking values within 
an open interval. Let 5* be a weakly trapped surface in H, i.e. both of its expansions are 
non-positive: 0(f) < 0, 0(„) < (in particular, S could be a marginally trapped surface). 
Then S must in fact be a marginally trapped surface and must coincide with one of the St- 
This shows that the foliation of H by the St must be unique. As a corollary, consider as in 
Fig. [9] a MTT generated by MTSs St which lie in spacelike surfaces St. This is a situation 
common in numerical relativity where one uses the Et for time evolution, and one locates 
MTSs on them. If we were to use a different foliation, and in particular a particular spatial 
hyper-surface E' which does not coincide with any of the Et. If the intersection E' p|H is not 
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Figure 12: This situation is ruled out by the uniqueness results of Ashtekar and Galloway. If there 
is a dynamical horizon %, then there cannot be another dynamical horizon T-L' lying completely in 
the past domain of dependence D^{T-L). 



one of the St , then it cannot be a marginally trapped surface. This is illustrated in Fig. [7| 
for the Vaidya spacetime. The intersection of the non-symmetric spatial hyper-surface with 
the spherically symmetric dynamical horizon is the red surface, and it is not a marginally (or 
weakly) trapped surface. Thus, if we had a different foliation Ej^ and we locate MTSs on these 
spatial hyper-surfaces, then we would end up with a different dynamical horizon H' . We note 
again the dramatic difference for an isolated horizon where every spherical cross-section is a 
marginally outer trapped surface. 

The question then arises, how different can H' be from H. A partial answer is provided 
by the next result [_13J : There are no weakly trapped surfaces contained in the past domain of 
dependence D~ (H) (apart from those which make up H itself)^ In particular, this result rules 
out a dynamical horizon contained in (J-L) — T-i. Once again, this result is illustrated by 
the Vaidya example discussed earlier in Fig. |4] The non-symmetric surface S in this figure is 
partly inside and partly outside the spherically symmetric dynamical horizon. More generally, 
weakly trapped surfaces lying partly inside D~ (H) are not ruled out. 

6 Outlook 

As discussed in the introduction, ever since the discover of the Schwarzschild solution al- 
most a century ago, research in black holes has led to seminal developments in theoretical, 

^The past domain of dependence D~ (J-L) is the set of spacetime points p such that every future causal curve 
from p intersects "H 
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observational and computational physics. This includes the singularity theorems, black hole 
thermodynamics, the uniqueness theorems, the cosmic censorship hypothesis, black hole en- 
tropy calculations in various approaches to quantum gravity, the various astrophysical phe- 
nomena involving black holes, and the recent results from numerical relativity. The framework 
of quasi-local horizons, which takes trapped and marginally trapped surfaces as its starting 
point, provides a unified approach for studying various aspects of black hole physics. While 
we have only touched upon a few aspects and applications of this framework, the material 
presented in this chapter will hopefully motivate the reader to delve further into the subject. 
An important theme in this discussion is that to base our understanding of black holes on 
lessons from stationary cases can be misleading. Dynamical situations have some essentially 
different features and intuition from stationary examples can easily lead us astray. We have 
illustrated this by the Schwarzschild and Vaidya examples. 

In this chapter we have introduced trapped surfaces and various kinds of quasi-local hori- 
zons through examples. The eventual goal of these studies (from a physics viewpoint) is to 
understand the properties of the surface of a black hole. We have discussed the inadequacy 
of event horizons for this purpose due to its teleological properties, and it is desirable to find 
a suitable replacement. Penrose's trapped surfaces and the boundary of the trapped region 
seem ideally suited for this task and lead naturally to the various definitions of quasi-local 
horizons. The simplest example is of course the Schwarzschild black hole. In this example, 
the boundary of the trapped region agrees with the event horizon, and both notions give rise 
to the same physical ideas. Difficulties arise however when we consider non-stationary black 
holes. We illustrated this through the imploding Vaidya spacetime. The intuitively obvious 
horizon, the analog of the r = 2M hyper-surface in Schwarzschild is a spherically symmetry 
dynamical horizon and it is separated from the event horizon. However, the non-spherically 
symmetric outer trapped surfaces extend up to the event horizon. To make matters more 
complicated, trapped surfaces do not extend all the way to the event horizon. We used these 
examples as motivations for general definitions and we reviewed some basic results regarding 
trapped surfaces and quasi-local horizons. We saw that marginally trapped surfaces are not 
as ill-behaved as one might think, and under physically reasonable conditions, they do evolve 
smoothly. The equilibrium case, described by isolated horizons is also of great interest. It 
covers a wide variety of situations where a black hole is in equilibrium in a dynamical space- 
time and is the best understood quasi-local horizon. Finally, in the dynamical case, we saw 
that one can assign physical quantities such as mass, angular momentum and fluxes through 
dynamical horizons. 

There are a number of topics that we have not discussed. In particular, we have not 
discussed the various applications of these notions in numerical relativity. Similarly, our dis- 
cussion has been restricted to the classical world and we have not talked about e.g. the 
quantization of isolated horizons and the black hole entropy calculations. Even in the dis- 
cussion of the mathematical properties of quasi-local horizons, we have discussed black hole 
multipole moments, symmetries, and inclusion of various kinds of matter fields only very 
briefly. Another significant omission is the discussion of black holes near equilibrium. Reviews 
of these topics can be found in e.g. [361 1241 116| 

We have seen that the outstanding question in this field is the non-uniqueness of dynamical 
horizons. We have discussed some restrictions on where dynamical horizons can be located. 
However, it is a fact that there are a multitude of smooth marginally trapped tubes and 
dynamical horizons in a general black hole spacetime, and there seems to be no obvious way 
of picking a preferred one. While there could be reasonable choices in specific examples, there 
does not seem to be any general solution to this problem. While one can get away with using 
event horizons in globally stationary spacetimes, choosing to live with event horizons in general 
is not an option because of its global properties. One could perhaps consider dealing with the 
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full set of dynamical horizons and marginally trapped tubes. Many of them are smooth one 
can study them individually; e.g. the dynamical horizon flux formulae apply to all dynamical 
horizons equally. However, if we wish to assign physical properties to the black hole such as 
mass, angular momentum, fluxes and higher multipoles, which one should we choose? We 
will generally get different results depending on our choice. The boundary of the trapped 
region could be a reasonable alternative, but as we have seen, it is generally not a marginally 
trapped tube itself. It is also not clear whether one can use this boundary to study, say, black 
hole thermodynamics and other physical phenomena that we believe are true for black holes 
and besides, this boundary also has a number of global properties and is difficult to locate 
(thus making it not better than event horizons in many regards). An interesting possibility, 
suggested by Bengtsson and Senovilla, is the core of the black hole region, i.e. the portion 
of the trapped region where all trapped surfaces penetrate. Removing the core would then 
completely eliminate all trapped surfaces. If it is indeed a proper subset of the trapped region, 
then is its boundary a dynamical horizon? There are indications that the region r < 2M{v) 
region of Vaidya is such a core, and its boundary is the spherically symmetric dynamical 
horizon. However the core is not unique and some cores are not spherically symmetric [23| . It 
is an interesting open question whether this idea can be developed further. 
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